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I. INTRODUCTION 

In previous work we presented theoretical arguments and strong phenomenological 
evidence that Regge trajectories of ordinary quark-antiquark mesons are essentially nonlinear 
and can be well approximated, for all practical purposes, by a specific function, the so 
called square-root form. With a few additional assumptions intended to reduce the number 
of independent parameters, and when possible, tested for self-consistency, we obtained a 
remarkable agreement with both bound (resonant) state (t > 0) and scattering (t < 0) data. 
What makes this success even more impressive is the fact that the input parameters, with 
the exception of one, were determined by the masses of only a few lowest lying bound states. 

The theoretical motivation for our previous study was the view that the properties of 
the gluon field, e.g. the flux tube, change with the increasing size of the hadron. At long 
enough distances, a linear potential is simply a wrong approximation to the interaction of 
meson constituents. Therefore, Regge trajectories cannot be asymptotically linear Q in t, 
even if they appear to be so over a limited range of t. In searching for a more practical 
approximation, we studied nonlinear Regge trajectories corresponding to dual amplitudes 
with Mandelstam analyticity (DAMA) 3[ in a simplified situation, that is, toy models 
which, nevertheless, maintain a .e^blJe to QCDflfla. Aft. though nation, 
we argued that (i) DAMA amplitudes can be expected to fit spectra; (ii) nonlinearity arises 
due to the color screening of the flux tube, and thus, (iii) one can expect the same qualitative 
behavior of the trajectories regardless of quantum numbers of the quark-antiquark meson. 

In this paper we attempt to go beyond ordinary quark-antiquark mesons to trajectories 
for glueballs. However, it is not a priori clear that non-gg mesons such as glueballs or hybrids 
can be satisfactorily described by the same form of Regge trajectories as are ordinary mesons. 
Even though the nonlinearity of Regge trajectories is due to color screening, which is not 
exclusive to qq mesons, the behavior of glue in exotic systems can be different. In particular, 
in the case of glueballs there have been contradictory opinions about the basic nature of 
these states, ranging from solitonic through loops of glue to glue strings. The latest lattice 
QCD results indicate that the loops-of-glue picture does not agree with the lattice spectra to 
the extent that the bag model does, thus supporting a constituent-like picture for gluons (in 
addition to quarks) 6]. (However, it should be kept in mind that the flux tube model is, by 
definition, a very simple model with a big symmetry group and so cannot be used to make 
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detailed predictions for the spectrum; the flux tube-like distribution of the fields should be 
measured instead.) The same authors view hybrids as states where a quark and an antiquark 
are subject to a potential corresponding, roughly speaking, to an excited flux tube. They 
generate a plethora of these potentials in quenched QCD. Obviously, in unquenched QCD 
the picture can be drastically different. It is, unfortunately, challenging and difficult at 
present to capture the physics of string breaking accurately with lattice methods f|. Here 
we focus on pure gluonic states and only comment on the prospect of extending our analysis 
to hybrid mesons in the conclusion. 

Encouraged by the lattice QCD analysis of Ref.|f|, we apply the same phenomenological 
approach that worked so remarkably well for ordinary mesons to the pomeron trajectory. In 
a departure from our previous calculations, we analyze scattering data to fit parameters of 
the square-root Regge trajectory. Naturally, in this case the extracted value of the threshold, 
or termination point of the real part of the trajectory, is subject to a large uncertainty. To 
reduce the uncertainty, we add to our data set the mass of the lowest lying tensor glueball 
determined in Ref.0. 

(There are a number of other lattice QCD determinations of this mass, see e.g. Refs.jlj], 
but the rms variation of the values reported is only about 1% which is considerably smaller 
than the uncertainties in the individual determinations. Hence, choosing one specific case 
has a negligible effect on our results compared to the overall (statistical and systematic) 
uncertainties below and avoids a perhaps unrealistic reduction of the total uncertainty.) 

With this addition, the uncertainties in the trajectory parameters are significantly reduced 
while the values are little affected. We take the value of the threshold as an indication of a 
maximum mass beyond which no glueball states exist, bearing in mind that we are unable 
to prove at present that the square-root form will be as comparably efficient in describing 
glueballs as it was for ordinary mesons. 

We should also note that the threshold we find here is larger than what is inferred from 
trajectories fitted to quenched lattice QCD glueball mass states. Whether this is due to 
the difference between quenched and unquenched QCD, or the functional forms used for the 
trajectories remains unknown. 

The choice of a specific trajectory within the allowed range of DAMA trajectories intro- 
duces an unknown systematic error. In case of ordinary mesons, our results justified the 
assumptions a posteriori. Due to insufficient data, this is not possible when dealing with 
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pure glue trajectories. To estimate the systematic errors, we repeat the fit with the other 
limiting form of DAMA trajectories, the so called logarithmic form. We take the difference 
between the thresholds obtained is indicative of the systematic errors. 

While the pomeron trajectory is the only glueball trajectory for which unambiguous 
data exist, it may be, unfortunately, affected by the gluon condensate. In the absence of 
experimental glue bound state data, we use glueball masses from lattice QCD to determine 
the thresholds of other, less peculiar trajectories, bearing in mind that the masses are subject 
to the quenched approximation. 

This paper is organized as follows: Section II starts with a brief introduction to the 
DAMA trajectories, followed by our fits to pomeron data. We compare the fits for the 
two limiting forms of trajectories, and discuss the physical meaning of the results. Section 
III is devoted to glueball spectroscopy. We conclude with a comment on possible future 
applications. 



II. DAMA TRAJECTORIES AND FITS TO DATA 



The class of dual models called dual amplitudes with Mandelstam analyticity (DAMA) 



is a generalization of Veneziano amplitudes [8| to the most general form consistent with 
Mandelstam analyticity. (DAMA has the Veneziano limit a(t) ~ t, but the transition to 
this limit occurs discontinuously |i|.) 

A meson trajectory (Xjj(t), can be parametrized on the entire physical sheet in the fol- 
lowing form: 

a fl {t)=a fl {Q)+ lu [T»--{T fl -ty], < v < ~. (1) 

(up to a power of a logarithm), assuming that: ctji(t) is an analytic function having a physical 
cut from some value t to oo; it is polynomially bounded on the entire physical sheet, and 
there exists a finite limit of the trajectory phase as \t\ — > oo The subscripts i, j indicate 
dependence of the parameters on the flavor content of the meson within a meson multiplet. 
In this paper we drop the subscripts for simplicity. 

The parameter 7^ is the universal asymptotic slope for nonlinear trajectories jll|. 

a(t) ~ -lv{-t) v , \t\ -> 00; 
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both 7„ and the exponent v are independent of quantum numbers. In order for the slope to 
be positive at small t, 7^ > 0. The intercept ctji(0) varies for different trajectories, and in 
accord with the Froissart bound should satisfy <x^(0) < 1. In reality, however, there is an 
exception — the intercept of the pomeron trajectory is observed to be slightly larger than 
1. The parameter T~j is often called the trajectory threshold. 
Note that for \t\ <T, Eq. (JTJ) reduces to a (quasi) linear form: 

a fi (t) = a fi {0) + vjT»fH = a fi {0) + a' fi {0) t. (2) 

The value of v is restricted to lie between and 1/2, in accordance with Ref. The 
value v = should be understood as a limit v — > 0, j^u fixed. In this limit, the difference 
of fractional powers reduces to a logarithm, viz., 

a(t) =a(0) -TioglogM- — J , 7log = lim 7 ^ (3) 

V J log/ 

Unlike a trajectory with any value of v 7^ 0, the real part of the "logarithmic" trajectory 
does not freeze-out when t reaches T. The real part continues to grow; the only change for 
t > T is that the trajectory acquires a constant imaginary part. 

The upper bound on v gives the so-called "square-root" trajectory, viz. 

a(t) = a(0) + 71/2 [VT-VT-t]. (4) 

When t reaches T, the real part of the "square-root" trajectory stops growing, and there 
are no states with a higher angular momentum than l m ax = &{T) ■ For this reason, the 
parameter T is also called the trajectory termination point. This is true for any value of 
v=£ 0. 



A. Fit with the square root form 

In Ref. jll we used only the square- root form of Eq.(j3} for spectroscopy purposes. We 
determined the value of 71/2 from the p trajectory, and then it was taken as universal for all 
other meson trajectories. Our calculation is in excellent agreement with various data, not 
only spectroscopic but scattering as well, and is self-consistent, justifying the assumptions 
a posteriori. This leads us to believe that the extracted value of 71/2 is reliable, and that 
the square-root form is close to the true functional form of meson trajectories. Since 7^ has 
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Pomeron Data and Square Root Fit 

with and without lattice glueball 
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- Fit to Pomeron data only 

- Fit with M(2++) glueball 
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FIG. 1: Pomeron data together with our best fits of the square-root form: The dashed line fit 
includes the averaged mass M(2 ++ ) as a data point; for the solid line fit only the scattering data 
was used. 
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to be a constant independent of the flavor content or quantum numbers of the meson, we 
apply the same value, 

7i /2 = 3.65 ±0.05 GeV" 1 , (5) 

to the data for the pomeron. 

With the value of the universal asymptotic slope fixed, we can now fit the pomeron 
scattering data to find the remaining parameter of the pomeron trajectory, T\/2- Of the two 



data sets |12j shown in Fig. 1, we use only 7p — > <pp since at ZEUS energies, it is devoid of 
significant contributions from exchanges other than the pomeron. For comparison, we show 
both sets of data. Note that the data points corresponding to ^p — > p°p are consistently 
above the jp —>■ <pp as expected due to contributions from additional exchanges. 
In Fig. 1, the solid line shows our best fit of the square-root form, Eq.(0J): 

a V2 (0) = 1.08 ±0.01, (6) 



Ti /2 = 11.56 ±2.08GeV (7) 

X 2 /d.o.f. = 6.07/9. 

The relatively large, 18% error on the fitted threshold is caused in part by the data 
uncertainties, and in part by the fact that all data points are concentrated in a small region 
of t near zero; consequently, \t\ << T. The fit is essentially dominated by the first term in 
the Taylor expansion of Eq. (j3J) : 

a 1/2 (t) = a (0) + -^=t (8) 



= 1.08 ±0.1578 [GeV" 2 ] t (9) 

In order to determine the threshold with a better accuracy, additional data is needed. 
We use the mass of the lowest lying tensor glueball from Ref. M(2 ++ ) = 2.40 ±0.13 GeV. 
With this additional data point, the error on the extracted threshold is reduced by factor of 
two, while the y 2 of the fit remains comparably small and the fitted parameter values are 
essentially unchanged. (See Figure 1). In this way we obtain: 

ai /2 (0) = 1.081 ±0.007, (10) 



T 1/2 = 11.57 il.lGeV (11) 

X 2 /d.o.f. = 6.07/10. 



We note that the x 2 /d.o.f. in both cases is smaller than one would expect on general 
grounds. This is similar to the case observed in Ref. 13] for other trajectories, presumably 
for similar reasons. 



B. Fit with the logarithmic form 

Even though the logarithmic trajectory, Eq.(|3J), itself is not realistic, since its real part 
grows without bound (in contrast to any other trajectory of our nonlinear form, Eq.(J2J), 
with v ^ 0), it is useful to study because the true trajectory can lie anywhere between the 
two limiting forms. In addition to this reason, comparison of the fits with the two limiting 
forms can, to some extent, illuminate the issue of systematics. There is undoubtedly a 
systematic error associated with the choice of DAMA trajectories as the class of trajectories 
within which the true trajectory lies, and this we cannot estimate. There is an additional 
systematic error arising from the specific choice of v — 1/2 within the model and from the 
way we determine the parameters of the trajectories. It is this second uncertainty that we 
address in this section. 

To fit the data with a logarithmic form, we first need to determine the value of 7i og . In 
complete analogy with our calculation utilizing the square root form, the value of 7i og is 
determined from the p trajectory. We find 

7i og = 8.00 ± 0.34 . (12) 

Note that 7i og is dimensionless. 

Our best fit of the logarithmic form, Eq.Q, to the scattering pomeron data, 

aiog(0) = 1.08 ±0.01, (13) 

sjr^g = 7.09±0.64GeV (14) 

xVd.o.f. = 6.07/9 

is shown in Fig. 2. Also in Fig. 2, we show the fit to the scattering data plus mass of the 
tensor glueball, yielding 

aiog(0) = 1.079 ±0.007, (15) 

^JtQ = 7.23±0.34GeV (16) 

xVd.o.f. = 6.12/10. 
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The addition of the bound state data point leads to a 2% increase in y T\ og and almost a 
factor of two reduction in its error. Recall that, for the square root form, the threshold 
remains almost the same (less than 0.1% change). 

Again, not surprisingly, the fit is dominated by the linear term 

a log (t) = a(0) + ^t (17) 

J- log 

= 1.08 + 0.1590 [GeV~ 2 ] t. (18) 



C. Estimate of systematic error 

The linear term dominance in all of our fits implies that the thresholds of the square root 
and logarithmic trajectories, extracted from pomeron data, are simply related. The relative 
size of the threshold for the two limiting forms follows directly from comparison of the two 
linearized trajectories, Eq. (JSJ) and Eq. (fT7j) : 

, 7iog 



T log = 2^JT 1/2 . (19) 

7l/2 

The ratio of is, in our calculation, fixed by the p trajectory data used as input. To find 
the three parameters of the p trajectory, (its intercept, threshold and the universal parameter 
7„ for any chosen v), we restrict the trajectory to pass through the three experimentally 
well established points, specifically the intercept, and the mass and spin of the p and the p$. 

The straight line that crosses p and p% leads to an intercept smaller that the observed 
value. This means that the linear form is insufficient Q|. By an a posteriori comparison 
of the DAMA trajectory to its truncated Taylor series, one can see that the fit is basically 
dominated by terms up to 0(t 2 ), viz. 

fi 1/2 (t) ~ «(0) + -^t + J^t 2 , (20) 

ZI l/2 ° I l/2 



a log (t) ~ a(0)+ (21) 
1 w 1 1 log 



(We use the tilde to distinguish the p trajectory from the pomeron trajectory discussed so 
far.) For example, the square root trajectory evaluated at the mass of the P3 differs from its 
Taylor series, Eq. (j20|) . by less than 0.5%. 

Setting the Taylor series coefficients in Eqs.(j2DJ) an d (|2*H) equal, we obtain 

Tio g ~ 2f 1/2 (22) 



7io g ^ 7i/2 (23) 

This leads the following relation for the pomeron thresholds: 

T log ~2 v^^T^. (24) 

The numerical values we find are in a good agreement with these relations. 

Is there anything deep about these relations? Note that Eqs. (j22|) and (J2*3|) are a direct 
consequence of the fit being dominated by up to quadratic terms in t and our requirement 
that the trajectory pass exactly through the three input points. Thus, Eq. (j23|) cannot hold 
for a universal 7 unless all thresholds are identical. Alternatively, the three point restriction 
is too strong and/or only one specific value of v can be correct. 



D. Understanding the systematic error in terms of toy models 

To understand this issue further, we turn to our toy models Within the framework 
of our generalized string model, it is possible to reconstruct the potential from a Regge 
trajectory 0). Earlier we found that potentials corresponding to the square-root and 
to the logarithmic trajectories, respectively, are very close, when normalised to the same 
asymptotic value. Furthermore, they are also very close to a potential found from a fit to 
lattice data jlj] that we used in another toy model, which consists, basically, of a leading 
order Born-Oppenheimer (LOBO) approximation for a system of a very heavy quark and 
antiquark. In that toy model, we found that the spectrum could be equally well fitted by 
both limiting forms of the Regge trajectories, with nearly the same thresholds 

The difference between the toy model study and the situation at hand is the following: 
In the toy model we fit the data with many points, in effect optimizing parameters of the 
underlying potentials so that they produce similar results for a large number of bound states. 
In the fit to real data, we have instead a very few lowest lying bound states, and solve for the 
parameters of trajectories. Implicitly, we demand that the underlying potentials produce the 
same results at those input points. Since the points are the lowest lying bound states, this 
corresponds to "aligning" the potentials in the region relevant for the lowest lying bound 
states, which can, and should be expected to, lead to different asymptotic values. Since the 
asymptotic value of the potential is directly related to the threshold of the Regge trajectory 
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(at least in the toy model), this translates into the thresholds for square root and logarithmic 
trajectories being different. 

The difference between the extracted thresholds is thus an indicator of systematic errors. 
Combining errors in quadrature, we conclude that the threshold for the pomeron trajectory 
is (9.4 ± 1.6) GeV. 

III. CONCLUSIONS ON GLUEBALL SPECTROSCOPY 

As is obvious from Eqs.(J7|) and (jll|) . the square root form of the trajectory with the 
parameters fitted to scattering data alone gives the same mass predictions as the fit to both 
the scattering data and the tensor glueball mass, but with larger errors. The fit to both 
the scattering data and the mass of the 2 ++ glueball from lattice may give quite precise 
predictions for higher excited states, providing the lattice value is close to the true mass of 
the glueball. Note that our method works very well even for higher excited states. Using 
the fit, Eq. ljllj) . we obtain the following predictions for excited glueball masses: M(4 ++ ) = 
4.21 ± 0.21 GeV, M(6++) = 5.41 ± 0.28 GeV, and we obtain M(2++) = 2.38 ± 0.12 GeV, 
with the same central value obtained from purely scattering pomeron data. 

Based on our calculation we conclude that the threshold for tensor glueballs can be 
expected in the region no lower than 7 to 8 GeV and no higher than 11 to 12 GeV. From 
the fitted values of the thresholds for various meson multiplets, we know that the thresholds 
for the same flavors, but different multiplets, vary by less than 20%. We expect the same to 
be true for glueballs. 

Note that glueball threshold is much larger than what we found for ordinary (qq) light 
mesons. This is related to the smaller (local in t, not asymptotic, of course) slope for the 
pomeron trajectory but also devolves from our DAMA approach. From this higher value, 
we infer that a significantly larger number of higher states are available for glueballs than 
for mesons. 

It is also interesting to contemplate what the maximal value of J = J max may be for the 
states allowed. We recall from our previous work jl| that the square root trajectory tends 
to overestimate the growth near the threshold; that is, at any given mass it tends to predict 
slightly high angular momentum for states approaching the threshold. This was expected 
from model studies and further confirmed by fits to light meson spectra where sufficient data 
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Pomeron Data and Logarithmic Fit 

with and without lattice glueball 
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- Fit to Pomeron data only 

- Fit with M(2++) glueball 
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FIG. 2: Pomeron data together with our best fits of the logarithmic form: The dashed line fit 
includes the averaged mass M(2 ++ ) as a data point; for the solid line fit only the scattering data 
was used. 

were available. For example, for the 02 trajectory the square root form allowed for a J = 8 
state, whereas we concluded that J = 6 should actually be the last state on the trajectory. 
The glueball trajectory allows for very large values of J max , possibly as high as 36. The 
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specific value of J max is obviously subject to a large uncertainty, but our conclusion is firm 
that the value will be significantly larger than that for ordinary mesons, possibly even larger 
than that for heavy quarkonia. This raises the paradoxical possibility that high J glueballs 
may be sufficiently narrow to identify experimentally. (The only other states expected to 
be available to mix with them would be heavy quarkonia. Such mixing is suppressed to 
perturbative values by the heavy quark mass.) 

Unfortunately, the pomeron trajectory is the only glueball trajectory for which experi- 
mental data is available, and we expect that there are certain exceptional aspects associated 
with it: As for the light quark-antiquark system, a perturbative analysis shows a strong 
attraction between two gluons with a threshold at zero. Relativistically, this suggests that 
the two particle bound state will develop a negative mass-squared, requiring 15j the for- 
mation of a gluonic vacuum condensate (which is known to occur) and a mass gap for the 
lowest scalar state above the shifted vacuum. The pomeron trajectory does pass through 
J = at a negative mass-squared and there is, of course, no physical state there. If the 
lowest mass scalar glueball is indeed the scalar state as shifted by the formation of the gluon 
condensate, it need not appear at the mass expected from standard consideration of the 
daughter trajectories of the pomeron. 

Such a distortion, however, does not appear to occur [l^ for the /o(980) trajectory in the 
case of light quarks, where the issues of 'four quark' states, as well as mixing (including with 
the scalar glueball), also arise. However, one does not know if this apparent regularity will 
hold for the pomeron trajectory states as well. Therefore, it is also interesting to investigate 
additional glueball states. 

In the absence of other data, we turn to lattice QCD. So far, only the spectrum of glueballs 
in quenched QCD has been calculated. Of the states listed in Ref.[6], the 0~ + at 2.59 ±0.17 
GeV and 2 _+ at 3.10 ± 0.18 GeV should form a common trajectory with an intercept 
—3.58 ± 1.48 and threshold y/T = 3.90 ± 0.91 GeV. Another trajectory is formed by a 1 H 
at 2.94±0.17 GeV and a 3 +_ at 3.55±0.21 GeV. Its intercept is -2.60±1.45 and threshold 
y/T = 4.87 ± 1.25 GeV. Note the large errors of the extracted values. The remaining data 
from Ref.Ja] are ambiguous for our purposes because they can contain admixtures of higher 
angular momentum states. 

The values of thresholds extracted from lattice data are significantly lower than the 
threshold of the pomeron trajectory, although, within the large errors, they are consistent 
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with the lower limit of the pomeron threshold. Moreover, since the thresholds for different 
trajectories of the same flavors can be expected to differ as much as 20%, the apparent 
disagreement is not alarming. Unquenched lattice data, therefore, could be very useful for 
further advancing our understanding the nature of the pomeron trajectory, in the absence 
of further experimental evidence. 



IV. FUTURE PROSPECTS 



As a prospect for future work, we comment on the possibility of applying our approach 
to hybrid mesons, which are of current interest to both experiment and theory. In Ref.0, 
the authors used the LOBO approximation to calculate spectra of bottomonium as well as 
bottomonium-like hybrids. The LOBO approximation was demonstrated to be an efficient 
and reliable method. However, the potentials used as inputs are generated in quenched 
lattice QCD, and thus, not all of the states predicted may survive in an unquenched theory. 
Color screening due to light quarks can be expected to become more and more important 
with increasing size of the hadron. 

Exactly how many of the hybrid states can be reliably extracted from a quenched calcu- 
lation is unclear. For example, the wave function of the lowest lying hybrid is found to be 
larger than that of the lowest lying quarkonium, but still smaller than the scale where flux 
tube breakage is expected There is evidence that survival of the lowest lying hybrids as 
well-defined resonances remains conceivable 

Unfortunately, implementing the physics of flux tube breaking in lattice QCD is very 
difficult at present. This is where our phenomenological approach can be of assistance. 
The key assumption is that because the curvature of the Regge trajectory arises due to the 
screening and breakage of flux tube, we can use the same (that is, square root) form of the 
trajectories as for ordinary mesons. 

It is also unfortunate that, at present only spin averaged lattice data are available for 

n 

the bottomonium-like hybrids [6]. If the physical states are nearly degenerate, i.e., spin 
splittings are small, then the spin averaged data can be used to extract the parameters of 
the Regge trajectories, and, in particular, the threshold value can be reliable. For example, 
in pure qq mesons, the bottom mass provides a sufficient suppresion factor for the subleading 
splittings, whereas the charm mass is insufficient. In case of bottomonium-like hybrids, most 
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possible spin dependent operators are suppressed by the heavy quark mass. However, there 
is a contribution due to the angular momentum of the glue that can be expected to be of the 
same magnitude as the LOBO splittings. This invalidates any conclusions one could draw 
using our approach based on the leading order splittings at present. Given more precise 
lattice data, it would be interesting to compare trajectory parameters of hybrids to those of 
ordinary mesons and or glueballs. 

It should also be worthwhile to examine daughter trajectories using our approach. Con- 
versely to the above, a plethora of relevant experimental evidence is available. A key question 
here is whether the thresholds for the daughter trajectories are consistent with those found 
for the parents. 

This research is supported by the Department of Energy under contracts W-7405-ENG-36 
and DE-FG02-87ER40365. 
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